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Abstract 



^ ' In the framework of the theory of open systems based on completely positive 

^ ■ quantum dynamical semigroups, we study the continuous variable entanglement 

O^i for a system consisting of two independent harmonic oscillators interacting with 

a general environment. We solve the Kossakowski-Lindblad master equation 
for the time evolution of the considered system and describe the entanglement 
in terms of the covariance matrix for an arbitrary Gaussian input state. Us- 
ff^ I ing Peres-Simon necessary and sufficient criterion for separability of two-mode 

' Gaussian states, we show that for certain values of diffusion and dissipation 

coefficients describing the environment, the state keeps for all times its initial 
type: separable or entangled. In other cases, entanglement generation, entan- 
Q>^ I glement sudden death or a periodic collapse and revival of entanglement take 

■ place. We analyze also the time evolution of the logarithmic negativity, which 

characterizes the degree of entanglement of the quantum state. 

X 

5—1 

5^ ; 1 Introduction 

The rapid development of the theory of quantum information has revived the interest 
in open quantum systems in connection, on one side, to the decoherence phenomenon 
and, on the other side, to their capacity of generating entanglement in multi-partite 
systems interacting with their environments. Quantum entanglement represents the 
physical resource in quantum information science which is indispensable for the de- 
scription and performance of such tasks like teleportation, superdense coding, quan- 
tum cryptography and quantum computation [1] . Therefore the generation, detection 
and manipulation of the entanglement continues to be presently a problem of intense 
investigation. 
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When two systems are immersed in an environment, then, besides and at the 
same time with the quantum decoherence, the environment can also generate a quan- 
tum entanglement of the two systems and therefore an additional mechanism to cor- 
relate them [21 [3] . In certain circumstances, the environment enhances entanglement 
and in others it suppresses the entanglement and the state describing the two sys- 
tems becomes separable. The structure and properties of the environment may be 
such that not only the two systems become entangled, but also such that a certain 
amount of entanglement survives in the asymptotic long-time regime. The reason is 
that even if not directly coupled, the two systems immersed in the same environment 
can interact through the environment itself and it depends on how strong this indirect 
interaction is with respect to the quantum decoherence, whether entanglement can be 
generated at the beginning of the evolution and, in the case of an affirmative answer, 
if it can be maintained for a definite time or it survives indefinitely in time [2]. 

In this work we study, in the framework of the theory of open quantum systems 
based on completely positive dynamical semigroups, the dynamics of the continuous 
variable entanglement for a subsystem composed of two identical harmonic oscillators 
interacting with its environment. We are interested in discussing the correlation 
effect of the environment, therefore we assume that the two systems are independent, 
i.e. they do not interact directly. The initial state of the subsystem is taken of 
Gaussian form and the evolution under the quantum dynamical semigroup assures 
the preservation in time of the Gaussian form of the state. 

The organizing of the paper is as follows. In Sect. 2 the notion of the quan- 
tum dynamical semigroup is defined using the concept of a completely positive map. 
Then we give the general form of the Kossakowski-Lindblad quantum mechanical 
master equation describing the evolution of open quantum systems in the Marko- 
vian approximation. We mention the role of complete positivity in connection with 
the quantum entanglement of systems interacting with an external environment. In 
Sec. 3 we write and solve the equations of motion in the Heisenberg picture for two 
independent harmonic oscillators interacting with a general environment. Then, by 
using the Peres-Simon necessary and sufficient condition for separability of two-mode 
Gaussian states [H [5] , we investigate in Sec. 4 the dynamics of entanglement for the 
considered subsystem. In particular, with the help of the asymptotic covariance ma- 
trix, we determine the behaviour of the entanglement in the limit of long times. We 
show that for certain classes of environments the initial state evolves asymptotically 
to an equilibrium state which is entangled, while for other values of the parameters 
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describing the environment, the entanglement is suppressed and the asymptotic state 
is separable. We analyze also the time evolution of the logarithmic negativity, which 
characterizes the degree of entanglement of the quantum state. A summary and 
conclusions are given in Sec. 5. 

2 Axiomatic theory of open quantum systems 

The time evolution of a closed physical system is given by a dynamical group Ut, 
uniquely determined by its generator H, which is the Hamiltonian operator of the 
system. The action of the dynamical group Ut on any density matrix p from the set 
Vil-i) of all density matrices in the Hilbert space Ti of the quantum system is defined 
by 

p{t) = Ut{p) = e-^^Ve^""* (1) 

for all t G (— cxD,c)o). According to von Neumann, density operators p G 'Dili.) are 
trace class (Tr p < oo), self-adjoint (p"'' = p), positive (p > 0) operators with Tr p = 1. 
All these properties are conserved by the time evolution defined by Ut- 

In the case of open quantum systems, the time evolution $t of the density opera- 
tor p[t) = $f(p) has to preserve the von Neumann conditions for all times. It follows 
that $4 must have the following properties: 

(i) <l>t(Aipi + A2P2) = Ai<l>t(pi) + A2$t(p2) for Ai, A2 > 0, Ai + A2 = 1, i. e. $j 
must preserve the convex structure of ViH), 

(ii) $t(pt) = <ft(p)^ 

(iii) $t(p) > 0, 

(iv) Tr <l>t(p) = 1. 

The time evolution Ut for closed systems must be a group Ut+s = UtUg. We have 
also f/o(p) = p and Ut{p) —>■ p in the trace norm when t — > 0. The dual group Ut 
acting on the observables A E B{T-C), i.e. on the bounded operators on H, is given by 

U^(A) = e^^*Ae-^^*. (2) 

Then Ut{AB) = Ut{A)Ut{B) and Ut{I) = I, where / is the identity operator on 
Ti. Also, Ut{A) — > A ultraweakly when t — > and Ut is an ultraweakly continuous 
mapping [SI El El- These mappings have a strong positivity property called complete 
positivity: 

Y,B}Ut{AlA,)B,>0, A,B,eB{n). (3) 
hi 
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In the axiomatic approach to the description of the evolution of open quantum 
systems [HllTllH], one supposes that the time evolution of open systems is not very 
different from the time evolution of closed systems. The simplest dynamics $j which 
introduces a preferred direction in time, characteristic for dissipative processes, is 
that in which the group condition is replaced by the semigroup condition [Gj [Qj [10] 

<^t+s = ^t^s, t,s>0. (4) 

The complete positivity condition has the form: 

Y^Bj^tiAlA,)B,>0, A,B,eBin), (5) 

where $t denotes the dual of $t acting on B{T-C) and is defined by the duality condition 

Tr($i(p)A) = Tr(p$i(A)). (6) 

Then the conditions Tr$((p) = 1 and = / are equivalent. Also the conditions 

$f(y4) —>■ A ultraweakly when t —>■ and ^ p in the trace norm when t —>■ 0, are 
equivalent. For the semigroups with these properties and with a more weak property 
of positivity than Eq. ([5]), namely 

A>0^ ^tiA) > 0, (7) 

it is well known that there exists a (generally unbounded) mapping L - the generator 
of $t, and $t is uniquely determined by L. The dual generator of the dual semigroup 
is denoted by L: 

Tr(L(p)A) = Tr(pZ(A)). (8) 

The evolution equations by which L and L determine uniquely $t and $t, respectively, 
are given in the Schrodinger and Heisenberg picture as 

dMp) 
dt 

and 

dMA) 
dt 

These equations replace in the case of open systems the von Neumann-Liouville equa- 
tions 



mtip)) (9) 



mm. (10) 
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and 

respectively. For applications, Eqs. (Q and (fTOl) are only useful if the detailed struc- 
ture of the generator L{L) is known and can be related to the concrete properties of 
the open systems described by such equations. For the class of dynamical semigroups 
which are completely positive and norm continuous, the generator L is bounded. In 
many applications the generator is unbounded. 

According to Lindblad [8], the following argument can be used to justify the 
complete positivity of $4: if the open system is extended in a trivial way to a larger 
system described in a Hilbert space Ti® K, with the time evolution defined by 

Wt{A ®B) = ®B, Ae B{H), B e B{iq, (13) 

then the positivity of the states of the compound system will be preserved by Wt 
only if $t is completely positive. With this observation a new equivalent definition 
of the complete positivity is obtained: is completely positive if Wt is positive for 
any finite dimensional Hilbert space K,. The physical meaning of complete positivity 
can mainly be understood in relation to the existence of entangled states, the typ- 
ical example being given by a vector state with a singlet-like structure that cannot 
be written as a tensor product of vector states. Positivity property guarantees the 
physical consistency of evolving states of single systems, while complete positivity 
prevents inconsistencies in entangled composite systems and therefore the existence 
of entangled states makes the request of complete positivity necessary 

A bounded mapping L : E{n) E{n) which satisfies Z(/) = 0, L{A^) = L(A)t 

and 

L{A^A) - L{A^)A - A^L{A) > (14) 
is called dissipative. The 2-positivity property of the completely positive mapping 

^t{A^A) > $i(At)$,(A), (15) 

with equality at t = 0, implies that L is dissipative. Conversely, the dissipativity 
of L implies that is 2-positive. L is called completely dissipative if all trivial 
extensions of L to a compound system described by ?i(S>/C with any finite dimensional 
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Hilbert space K, are dissipative. There exists a one-to-one correspondence between 
the completely positive norm continuous semigroups $t and completely dissipative 
generators L. The following structural theorem gives the most general form of a 
completely dissipative mapping L [8]. 

Theorem. L is completely dissipative and ultraweakly continuous if and only if 
it is of the form 

L{A) = '-[H,A] + ^E(^'[A^.] + [y^m), (16) 

where Vj, V^^V, e B{n), H E B{n)^ ^ . 

The dual generator on the state space (Schrodinger picture) is of the form 

L{p) = -^[H,p] + l.J2{[V^p,vl] + [V„pVj]). (17) 
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Eqs. ([9]) and (fTTIl give the explicit form of the Kossakowski-Lindblad master equation, 
which is the most general time-homogeneous quantum mechanical Markovian master 
equation with a bounded Liouville operator [HI [IDl [HI [12]: 

= -^[H, Mp)] + ^ T^ilvMp)^ vj] + [v„ Mp)vJ]). (18) 

j 

The assumption of a semigroup dynamics is only applicable in the limit of weak 
coupling of the subsystem with its environment, i.e. for long relaxation times [13] . 
We mention that the majority of Markovian master equations found in the literature 
are of this form after some rearrangement of terms, even for unbounded generators. 
It is also an empirical fact for many physically interesting situations that the time 
evolutions $t drive the system towards a unique final state p{oo) = limj^oo ^t{p{0)) 
for all p(0) G V{n). 

3 Time evolution of two independent harmonic os- 
cillators interacting with the environment 

We are interested in the dynamics of entanglement in a subsystem composed of two 
identical non-interacting (independent) harmonic oscillators in weak interaction with 
a general environment, so that their reduced time evolution can be described by a 
Markovian, completely positive quantum dynamical semigroup. If is the dynamical 
semigroup describing the irreversible time evolution of the open quantum system 
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in the Heisenberg picture, then the Kossakowski-Lindblad master equation has the 
following form for an operator A (see Eqs. (fTOl). (fT6l)) P [TOl [TTl [12] : 



dt 



(19) 



Here, H denotes the Hamiltonian of the open system and the operators Vj,Vj , defined 
on the Hilbert space of H, represent the interaction of the open system with the 
environment. We are interested in the set of Gaussian states, therefore we introduce 
such quantum dynamical semigroups that preserve this set. Consequently H is taken 
to be a polynomial of second degree in the coordinates x,y and momenta Px,Py of 
the two quantum oscillators and Vj, are taken polynomials of first degree in these 
canonical observables. Then in the linear space spanned by the coordinates and 
momenta there exist only four linearly independent operators Vj=i_2,3,4 



(20) 



where a^j, dyj^b^j^byj E C. The Hamiltonian H of the two uncoupled identical har- 
monic oscillators of mass m and frequency lj is given by 

1 



H 



mio 



2m 



iPx + pt) + ^i^ +y) 



(21) 



The fact that is a dynamical semigroup implies the positivity of the following 
matrix formed by the scalar products of the four vectors a^, b^^, a^, hy, whose entries 
are the components axj,bxj,ayj,byj, respectively: 



/ (a^.aa;) {a^h^) (a^a^^) (a^b^) \ 



-h 



^y^x) 

V (bj,a^) {hyh 



Ayhx) 



(22) 



(b^-a^) {Khyj 

(hyay) (hyhy) ) 

Its matrix elements have to be chosen appropriately to suit various physical models of 
the environment. For a quite general environment able to induce noise and damping 
effects, we take this matrix of the following form, where all the coefficients D^x-, Dxp^,--- 
and A are real quantities, representing the diffusion coefficients and, respectively, the 
dissipation constant: 



/ Dxx 

-Dxp^ + thX/2 



-D 



'xy 



XPx 

-D 



ihX/2 



V 



-D 



XPy 



D 



ypx 

PxVy 



D. 
-D 
D 



xy 

ypx 



-D 



XPy 



yy 



-D 



-D 



ypy 



ih\/2 



PxPy 

yp^ - ihX/2 

DpyPy 



(23) 
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It follows that the principal minors of this matrix are positive or zero. From the 
Cauchy-Schwarz inequality the following relations for the coefficients defined in Eq. 
( !23|) hold (from now on we put, for simplicity, h = 1): 



DxxDp^p^ 

DxxDyy 



X2 

>- 

XPx — ^ 



^yy^PvPv 



r.2 

D > — 

ypy — 4 ' 



DXXDpyPy 



D% > 0, Dp^p^Dp^py - Dl^^ > 0, 
> 0, DyyDp^p^ - > 0. 



XPy 



(24) 



The matrix of the coefficients fl23l) can be conveniently written as (T denotes the 
transposed matrix) 

in terms of 2 x 2 matrices Ci = Ci\ C2 = C-^ and C3. This decomposition has 
a direct physical interpretation: the elements containing the diagonal contributions 
C\ and C2 represent diffusion and dissipation coefficients corresponding to the first, 
respectively the second, system in absence of the other, while the elements in C3 
represent environment generated couplings between the two oscillators, taken initially 
independent. 

We introduce the following 4x4 bimodal covariance matrix: 



/ (^xxif) (^xp^{t) CFxy{t) 



a{t) 



a. 



^xpAt) 
it) 



a, 



PxPx 



xy 



(t) ayp^{t) 



XPy 



a, 



PxPy 



yy\ 
ypy^ 



cr„„,., (t) 



ypy 

PyPv 



it)\ 
it) 
it) 



(26) 



CTr,.,.,,., it) / 



ypx 

\(^XPyit) CTp^Pyit) 

with the correlations of operators Ai and A2, defined by using the density operator 
p of the initial state of the quantum system, as follows: 



^A.A^it) = -Tt{p{A,A2 + A2A,){t)) - Tr(pA(t))Tr(pA2(t)). 



(27) 



By using Eq. ( fT9i) we obtain by direct calculation the following systems of equations 
for the quantum correlations of the canonical observables [E 

da{t) 



dt 



Ya{t)+a{t)Y' +2D, 



(28) 



where 



Y 



-X 






1/m 
-A 












-A 



—muo 





1/m 
-A / 



(29) 
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D 





^xpx 


Dxy 


Dxpy 


\ 




^PxPx 


^ypx 


^PxPy 




Dxy 


^ypx 


^yy 


^ypy 




V ^xpy 


^PxPy 


^ypy 


^PyPy 


) 



(30) 



Introducing the notation a{oo) = limt^oo<^{t), the time-dependent solution of Eq. 
(I2HD is given by pi] 

a{t) = M{t){a{0) - a{oo))M^{t) + a(oo), (31) 

where the matrix M(t) = exp{Yt) has to fulfill the condition limf_^oo M{t) = 0. In 
order that this limit exists, Y must only have eigenvalues with negative real parts. 
The values at infinity are obtained from the equation 

Ya{oo) + a{oo)Y^ = -2D. (32) 

4 Dynamics of entanglement 

The two-mode Gaussian state is entirely specified by its covariance matrix (1261) . which 
is a real, symmetric and positive matrix with the following block structure: 



^(t) = u:^ " , (33) 




where A, B and C are 2x2 matrices. Their entries are correlations of the canonical 
operators x, y,Px and Py\ A and B denote the symmetric covariance matrices for the 
individual reduced one-mode states, while the matrix C contains the cross-correlations 
between modes. The elements of the covariance matrix depend on Y and D and can 
be calculated from Eqs. ([SID, (123) • Since the two oscillators are identical, it is 
natural to consider environments for which the two diagonal submatrices in Eq. fl2^ 
are equal, Ci = C2, and the matrix C3 is symmetric, so that in the following we 
take Dxx = Dyy, Dxp., = Dyp^, Dp^p^ = Dp^p^, Dxp^ = Dyp^. Then both unimodal 
covariance matrices are equal, A = B, and the entanglement matrix C is symmetric. 

4.1 Time evolution of entanglement 

It is interesting that the general theory of open quantum systems allows couplings 
via the environment between uncoupled oscillators. According to the definitions of 
the environment parameters, the diffusion coefficients above can take non-zero values 
and therefore can simulate an interaction between the uncoupled oscillators. Conse- 
quently, the cross-correlations between modes can have non-zero values. In this case 
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the Gaussian states with det C > are separable states, but for det C < it may be 
possible that the states are entangled, as will be shown next. 

In order to investigate whether an external environment can actually entangle 
the two independent systems, we use the partial transposition criterion [H [5] : a state 
is entangled if and only if the operation of partial transposition does not preserve its 
positivity. For the particular case of Gaussian states, Simon [5] obtained the following 
necessary and sufficient criterion for separability: S* > 0, where 

S = detAdetB + {^-\ det C\f - Tt[AJCJBJC^ J] - ^(det A + det B) (34) 

and J is the 2x2 symplectic matrix 

In the following we consider such environment diffusion coefficients, for which 

m'^uj^D^^ = Dp^p^, D^p^ = 0, m'^uj^D^y = Dp^p^. (36) 



This corresponds to the case when the asymptotic state is a Gibbs state 

In order to describe the dynamics of entanglement, we have to analyze the time 
evolution of the Simon function S(t) flM|) . We consider two cases, according to the 
type of the initial Gaussian state: separable or entangled. 

1) To illustrate a possible generation of the entanglement, we represent in Figure 
1 the function S{t) versus time t and diffusion coefficient D^py = d for a separable 
initial Gaussian state with initial correlations axxiS^) = 1, o"p^p^(0) = 1/2, crxp^{0) = 
0, (Jxy{0) = 0, o"p^py(0) = 0, (Txpy{0) = 0. We notice that, according to Peres-Simon 
criterion, for relatively small values of the coefficient d, the initial separable state 
remains separable for all times. For larger values of d, at some finite moment of time, 
when S{t) becomes negative, the state becomes entangled. In some cases the entangle- 
ment is only temporarily generated, that is the state becomes again separable. After 
that, at a certain moment of time, one can notice again a revival of entanglement. 
In these cases the generated entangled state remains entangled forever, including the 
asymptotic final state. 

2) The evolution of an entangled initial state is illustrated in Figure 2, where we 
represent the function S(t) versus time and diffusion coefficient d for an initial entan- 
gled Gaussian state with initial correlations (y^xiO) = 1, crp^p^(O) = 1/2, cTxp^iO) = 
0, CTxyi^) = 1/2, o"p^py(0) = —1/2, (TxpyiO) = 0. We notice that for relatively small 
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Figure 1: Simon separability function S versus time t and environment coefficient d, 
for A = 0.2, D = 0.115 and a separable initial Gaussian state with initial correlations 
a,,(0) = 1, a,^pM = 1/2, ^.pM = cT^yiO) = cr^^pyiO) = ^,.^^(0) = 0. We have 
taken m = uj = h = 1. 

values of d, at some finite moment of time S{t) takes non-negative values and there- 
fore the state becomes separable. This is the so-called phenomenon of entanglement 
sudden death. Depending on the values of the coefficient d, it is also possible to have 
a repeated collapse and revival of the entanglement. One can also show that for rel- 
atively large values of the coefficients D^^ and d, the initial entangled state remains 
entangled for all times. 

4.2 Asymptotic entanglement 

On general grounds, one expects that the effects of decoherence, counteracting en- 
tanglement production, is dominant in the long-time regime, so that no quantum 
correlation (entanglement) is expected to be left at infinity. Nevertheless, there are 
situations in which the environment allows the existence of entangled asymptotic 
equilibrium states. From Eq. (l32l) we obtain the following elements of the asymptotic 
entanglement matrix C(oo): 

""^y^"^)- 2m2A(A2 + c.2) ' ^^^^ 
cTxpyioo) = aypAoc) = 2m(A^^c.^) ' ^^^^ 
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Figure 2: Same as in Fig. 1, for an entangled initial Gaussian state with initial 
correlations a^^{0) = 1, (Tp^p,(0) = 1/2, cT^p,(0) = 0, cr^j^y^O) = 1/2, crp^py{0) = 
-1/2, a,p^(0) = 0. 

. , _ m^iJ^D.,y - 2miJ^XD,p^ + (2A^ + u;^)Dp^p, 

The elements of matrices A{oo) and B{oo) are obtained by putting x = in the 
previous expressions. We calculate the determinant of the entanglement matrix and 
obtain: 

det C(oo) 

= 4;^2(p\^2) X [i^^'D^y + ^Dp^vyf + 4A^(D.,Z),.,, - DlJ]. (40) 

With the chosen coefficients (!36i) . the Simon expression (IMl) takes the following form 
in the limit of large times: 

S{oo)=( ^^f ^+ ^% -- -4 41 

For environments characterized by such coefficients that the expression S{oo) fHT]) is 
strictly negative, the asymptotic final state is entangled. Just to give an example, 
without altering the general features of the system, we consider the particular case of 
Dxy = 0. Then we obtain that 5(00) < 0, i.e. the asymptotic final state is entangled, 
for the following range of values of the coefficient Dxpy characterizing the environment 
[I51[T6]: 

muDxx _ 1 Dxpy mujDxx 1 
A 2^^ ./WT^ A +2' 
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Figure 3: Logarithmic negativity L versus time t and environment coefficient d, for 
A = 0.2, D = 0.115 and a separable initial Gaussian state with initial correlations 



cr^^(O) = 1, ap^p,(0) = 1/2, (T^p,(0) = (r^y{0) = ap^p,(0) = (T^p^(O) = 0. We have 



taken m = uj = h = 1. 

where the diffusion coefficient D^^ satisfies the condition mujD^^/X > 1/2, equivalent 
with the unimodal uncertainty relation. We remind that, according to inequalities 
( !24|) . the coefficients have to fulfill also the constraint D^x > D^py- If the coefficients 
do not fulffi the inequalities then S'(oo) > and the asymptotic state of the 

considered system is separable. These results show that, irrespective of the initial 
conditions, we can obtain either an separable or an inseparable asymptotic entangled 
state, for a suitable choice of the diffusion and dissipation coefficients. 

4.3 Logarithmic negativity 

We apply the measure of entanglement based on negative eigenvalues of the partial 
transpose of the subsystem density matrix. For a Gaussian density operator, the 
negativity is completely defined by the symplectic spectrum of the partial transpose 
of the covariance matrix. The logarithmic negativity L{t) = — | log2[4/(cr(t))] de- 
termines the strength of entanglement for L{t) > 0. If L{t) < 0, then the state is 
separable. Here 




f{a{t)) = -(det A + det B) - det C 

.2 \i/2 
A + det B) - det C - det a{t) 




(43) 



13 



Figure 4: Same as in Fig. 3, for an entangled initial Gaussian state with initial 
correlations a^^{0) = 1, (Tp^p,(0) = 1/2, cT^p,(0) = 0, cr^j^y^O) = 1/2, crp^py{0) = 
-1/2, a,p^(0) = 0. 

In Figures 3 and 4 we represent the logarithmic negativity L{t) versus time t and 
diffusion coefficient D^p^ = d for the two types of the initial Gaussian state, separa- 
ble or entangled, with the same initial correlations, previously considered when we 
analyzed the time evolution of the Simon function S{t). As expected, we remark that 
the logarithmic negativity has a behaviour similar to that one of the Simon function 
in what concerns the characteristics of the state of being separable or entangled. De- 
pending on the values of the environment coefficients, the initial state can preserve 
for all times its initial property - separable or entangled, but we can also notice the 
generation of entanglement or the collapse of entanglement (entanglement sudden 
death) at those finite moments of time when the logarithmic negativity L{t) reaches 
zero value. One can also observe a repeated collapse and revival of the entanglement. 
In the case of an entangled initial state, the logarithmic negativity is a fluctuating 
function and decreases asymptotically in time. 

In our case the asymptotic logarithmic negativity has the form 



It depends only on the diffusion and dissipation coefficients characterizing the envi- 
ronment and does not depend on the initial Gaussian state. One can easily see that 
the double inequality (H2l) . determining the existence of asymptotic entangled states 
(S'(oo) < 0) is equivalent with the condition of positivity of the expression ( H^ of the 
logarithmic negativity, L{oo) > 0. 




(44) 
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5 Summary 



We have given a brief review of the theory of open quantum systems based on com- 
pletely positive quantum dynamical semigroups and mentioned the necessity of the 
complete positivity for the existence of entangled states of systems interacting with an 
external environment. In the framework of this theory we investigated the dynamics 
of the quantum entanglement for a subsystem composed of two uncoupled identical 
harmonic oscillators interacting with a common environment. 

By using the Peres-Simon necessary and sufficient condition for separability of 
two-mode Gaussian states, we have described the generation and evolution of en- 
tanglement in terms of the covariance matrix for an arbitrary Gaussian input state. 
For some values of diffusion and dissipation coefficients describing the environment, 
the state keeps for all times its initial type: separable or entangled. In other cases, 
entanglement generation or entanglement suppression (entanglement sudden death) 
take place or even one can notice a repeated collapse and revival of entanglement. We 
have also shown that, independent of the type of the initial state, for certain classes 
of environments the initial state evolves asymptotically to an equilibrium state which 
is entangled, while for other values of the coefficients describing the environment, the 
asymptotic state is separable. We described also the time evolution of the logarithmic 
negativity, which characterizes the degree of entanglement of the quantum state. 

The existence of quantum correlations between the two considered harmonic 
oscillators interacting with a common environment is the result of the competition 
between entanglement and quantum decoherence. From the formal point of view, the 
generation of entanglement or its suppression (entanglement sudden death) corre- 
spond to the finite time vanishing of the Simon separability function or, respectively, 
of the logarithmic negativity. Presently there is a large debate relative to the physi- 
cal interpretation existing behind these fascinating phenomena. Due to the increased 
interest manifested towards the continuous variables approach to quantum infor- 
mation theory, these results, in particular the possibility of maintaining a bipartite 
entanglement in a diffusive-dissipative environment for asymptotic long times, might 
be useful in controlling the entanglement in open systems and also for applications 
in the field of quantum information processing and communication. 
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